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Abstract 

The problem of an open string in background -B-field is discussed. Using the 
discretized model in details we show that the system is influenced by infinite number 
of second class constraints. We interpret the allowed Fourier modes as the coordi- 
nates of the reduced phase space. This enables us to compute the Dirac brackets 
more easily. We prove that the coordinates of the string are non-commutative at the 
boundaries. We argue that in order to find the Dirac bracket or commutator algebra 
of the physical variables, one should not expand the fields in terms of the solutions 
of the equations of motion. Instead, one should impose the set of constraints in 
suitable coordinates. 
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1 Introduction 



The constrained systems were first introduced by Dirac [1] within the discussion about 
singular Lagrangians. In such systems, Euler-Lagrange equations of motion lead to some 
acceleration- free identities as consequences of the singularity of Hessian (second derivative 
of Lagrangian with respect to the velocities). At the Hamiltonian level, this leads to 
emergence of constraints, i.e. functions of phase space coordinates which should vanish 
on- shell. 

We call that the first class constraints arc responsible for gauge transformations; while 
second class constraints restrict the system to a smaller sub-manifold of the phase space 
in which a Poisson structure can be recognized. In a very simplified picture, the first 
class constraints may be visualized as some momenta, so that they are involutive and 
generate transformations in their conjugate coordinates (i.e. the gauged variables). On 
the other hand, the second class constraints may be visualized as conjugate pairs of 
coordinates and momenta whose Poisson brackets in the original phase space arc nonzero. 
This seems to contradict the fact that they vanish on the surface of motion. Hence, it 
is necessary to introduce a new bracket, i.e. the Dirac bracket, such that the physical 
quantities (quantities defined on the constraint surface) have vanishing brackets with the 
constraints. A complete review on constrained systems can be found in [2] 

The singularity of the Lagrangian is not necessarily the only origin of the constraints. 
One may find constraints in the Lagrangian formalism in any way that one can impose 
acceleration-free equations on the system; or in the Hamiltonian formalism, in any way 
that one can impose a primary constraint on the system. As far as we know, apart from 
artificial problems in which primary Hamiltonian constraints are imposed by hand, no 
serious model has been introduced in which constraints emerge naturally, but not as the 
result of the singularity of the system. 

In the latest years of the last millennium, in another branch of research, a new phe- 
nomenon was discovered which has relationships with the constrained systems. It was 
found that for an open string coupled to a background 5-field, the canonical quantization 
procedure fails at the end points [3, 4, 5]. It was also observed in Refs. [6, 7] based on 
the work of Ref. [8], that in the presence of mixed boundary conditions due to i?-field, 
the propagators of coordinate fields possess singularities which can be interpreted as the 
non-commutativity of coordinate fields specially at the boundaries (branes). The origin of 
this non-commutativity is intensively discussed from then on [9, 10, 11, 12, 13, 14, 15, 16]. 

This observation led some authors [4, 5, 9] to deduce that the Dirac mechanism of 
second class constraints has some role at the boundaries of the string, hence the idea of 
considering the boundary conditions as Dirac constraints was born [4, 5] . Since the bound- 
ary conditions just put limitations on the solutions of equations of motion; if accepted as 
constraints, they should somehow be related to the second class systems. 

In the next section we will review the main set up of the above model, following mainly 
the method of reference [9]. Assuming the boundary conditions as constraints, we use the 



total Hamiltonian to impose their consistency. As we will see, the Lagrange multipliers 
are determined at the first stage of consistency (to be zero), but at the same time the 
second level constraints emerge. Then third level constraints are derived by demanding 
the consistency of second level ones, and so forth. In this way the procedure of finding 
the constraints proceeds unlimitedly. This suggests a new category in comparison with 
the ordinary theory of constrained systems in which the constraint chains terminate after 
a limited number of steps either by arriving at an identity (for first class systems), or by 
determining the Lagrange multipliers (for second class systems) [17]. 

During recent years several different and even opposite approaches have been applied 
to the problem, thoTigh most of them have derived similar results. Some authors do not 
consider the consistency of the constraints completely. They calculate the Dirac brackets 
just by using primary or at most second level constraints [10, 18, 19]. The idea of infinite 
number of second class constraints at the boundaries, however, was accepted by many 
authors [4, 5, 9], although fewer people note what happens to the Lagrange multipliers. 
Some authors have also tried to overcome the problem by considering discrctizcd model 
[9] or by following the symplectic approach [11, 20, 21]. Such approaches, though verify 
that the system is constrained in some sense, do not give final justification about the 
constraint characteristics of the system. 

Our first objective in this work is to find a distinctive understanding about the con- 
straint structure of the model. We want to know in what sense we find infinite number 
of constraints. We think that the key to study this part of the problem is the discretized 
model. It is known that the primary constraints emerge as the continuum limit of the 
equations of motion of the end points in the discretized model [9, 10]. However investi- 
gating the consistency conditions to deduce the set of secondary constraints has not been 
studied in discretized model yet. In sections (3) and (4) we have done this. We will see 
that considering the continuum limit together with the physical condition of continuity 
gives the set of desired constraints. We think that in the light of this study of discrete 
model any doubt about the existence of infinite number of constraints disappears. 

Next problem concerns the properties of the reduced phase space. Specially, one needs 
to know about the induced brackets on this space, which is the same as the Dirac brackets. 
With infinite number of constraints, the matrix of Poisson brackets of constraints is infinite 
dimensional. One needs to invert it to find the Dirac brackets. Some authors, which have 
accepted the existence of infinite number of constraints [5, 9], have tried to solve this 
problem directly. However, as we will explain in the following section, the mathematical 
manipulations given in different papers are not convincing yet. In fact, since people 
excepted to end up with non-commutativity in brane coordinates upon physical intuitions, 
most of the authors (except a few of them) have derived the famous result of Ref. [6] with 
more or less problematic mathematical methods. 

In section (5) we will try to give a very simple approach to find the Dirac brackets on 
the basis of Fourier modes. We show that this powerful physical tool, (considered in some 
different sense in Ref. [22]) serves as the suitable coordinates for describing the reduced 



phase space. Prom this point of view, after imposing the constraints, the brackets of the 
remaining physical modes can be written automatically; and then the Dirac brackets of 
the original fields can be derived, using their expansion in terms of physical modes. Apart 
from convenience in calculations, this method helps one understand better the mechanism 
which leads to unusual brackets of coordinate fields at the boundaries. 

In section (6) we give a comprehensive discussion about using the classical equations 
of motion of the fields in the process of quantization. We argue that in fact, the algebra 
of the observables is the essential entity in the quantization process rather than their 
dynamics. Wc show that just the dynamics of constraints is essential to be investigated 
before quantization. In this sense, we have derived the quantum properties of the reduced 
phase space, including the non-commutativity of brane coordinates, without using the 
solutions of the equations of motion. 



2 Problem setup 



Consider an open string with coordinate fields living in a target space specified by 
II = 0,1,.. .,D. Suppose that the string is coupled to a given antisymmetric field JF^,^ 
which for simplicity we assume no dynamics for it. The end points cr = 0, / of the string 
are constrained to move on a p + 1 dimensional Dp-brane characterized by X" = for 
a = p + 1, D. The U{1) gauge fields, Ai i = 0, 1, ...,p are also coupled to the string on 
the boundary (Dp-brane). The action of the model can now be written as 



1 



dadr 



G.,daX>'d,X''g-' + e'^'B.^daX'^d.X 



ab ; 



1 



drAidrX' 



(2.1) 



where 812 is the boundary of target space E, gab is the metric of world-sheet and e"* is the 
antisymmetric tensor on the world-sheet. We considered for simplicity only the bosonic 
sector. Similar arguments can be applied to a superstring [13, 23]. The bulk and boundary 
fields can be combined to construct modified Born-Infeld field strength = B — dA. Let 
also both end points attach to the same brane. Assume Oi' = ^ and the background 
metric to be fiat: 6*^,^ = Suppose, moreover, that the field strength is constant 

everywhere, given by the constant antisymmetric matrix Bij. Using these simplifications 
the action reads: 



dadr 



daX^'dbX^^g"' + BijdaX'dbXh 



ab 



fixing the diffeomorphisms and scaling invariance, the above action takes the form 



S 



dadr 



JCX^ - X'^X'„ + 2BijX'X'^ 



(2.2) 



(2.3) 



where "dot" and "prime" represent differentiating with respect to r and a respectively. 
Besides the conformal symmetry, this action has a global symmetry under transforming 
the coordinate fields by a constant, i.e. 



(2.4) 



Varying the action with respect to the fields gives the equations of motion 



^^X^'{a,T) -dlX''{a,T) ^0 // = 0,1,..,L> (2.5) 
together with the boundary conditions 

d,X'{a, r) + BijdrX^ia, r) = i = 0, 1, • • • ,p 

at the end points a = 0,1. As is apparent the boundary conditions in the directions 
perpendicular to the Dp-brane are simply of Dirichlet type; while along the Dp-brane we 
have mixed boundary conditions. This is the only place where the effect of the 5-field 
is experienced. Let concentrate in the following on this part of the problem and ignore 
the coordinate fields with Dirichlet boundary conditions. One simple way to realize this 
point is to assume that p — D. On the other hand, at the end of calculations one can 
turn off the S-field in as many direction as one desires to achieve results concerning the 
Neumann boundary conditions. The equations (2.6) do not contain accelerations, so it 
may be viewed as primary Lagrangian constraints. The canonical momentum fields and 
the Hamiltonian of the system read 



2 Jo 



Ui{a, r) = drX{a, r) + Bijd,Xj{a, r) (2.7) 

f UUi - Bijd^Xjf + {d^Xiy] da. (2.8) 
Jo L ^ 

The primary constraints (2.6) in terms of the phase space variables are 

= $,(a)|,=o, ^i = M<^)\a=i (2.9) 

with 

^i{a, t) = Mijd^Xjia, r) + B^jUjia, r), M^l-B^. (2.10) 
As for every constrained system, given the primary constraints (2.9) the total Hamiltonian 

Ht^H + \^] + ~X,^] (2.11) 

is responsible for the dynamics of the system. Hence, the consistency of the constraints 
gives 



$J = {$1, Ht} = H} + A,{$i, $]} = 
i>' = {$1, Ht} = H} + A,{<|i, $]} = 0. 



Using Eqs. (2.8-2.10) and the fundamental Poisson brackets 



{X,{a,T),X,{a',T)} = 

{X,{a, t), n,(a', r)} = 6,,6{a - a') (2.13) 
{n,(a,r),H,(a',r)} = 0, 



it is easy to see that 

^^^={<!?lH}^dMi\a=o (2.14) 
{$1, $]} = -2{MB)ij j 6{a)6{a')dJ{a - a')dada' . (2.15) 

Similar equations also arise for $j at the end point a = I. Let see how can the consistency 
conditions (2.12) come true. Noticing (2.14) and (2.15), wc find that the two Poisson 
brackets appearing on the right hand side of equations (2.12) are not of the same order. 
This fact is explained in more detail in [9] using the regularized form of the Dirac delta 
functions. Therefore the only way to satisfy the consistency conditions (2.12) , is to 
assume that 

A, = A, = 0, , , 

^ ! ' 2.16 

= $2 = 0. 

The important result is that the secondary constraints emerge while the Lagrange mul- 
tipliers arc determined. One should continue the consistency process by demanding the 
time derivatives of $^ and to vanish. From (2.11) and (2.16) we have Ht = H from 
now on, so 

$f = $2 = {^j^H} = a2$(a,r)U=o = (2.17) 
with similar expressions for $f. In this way two infinite constraint chains appear as 

$"=<( ^""=' (2.19) 

Next the problem arises: what are the Dirac brackets of the fields due to the above 
infinite constraints. We recall that the Dirac bracket of two quantities / and g in phase 
space is defined as 

{/, 9}d.b = {/, 9} + {/, X^}C'Hxj: 9} (2.20) 
where Xi ^-^e second class constraints and C^^ is the inverse of 

Cij^{Xi,Xj}- (2.21) 

In the present problem Cij is infinite dimensional and it is difficult (or in fact impossible) 
to find its inverse. It is obvious that for all m , n , i and j 

|$n^m|^Q (2.22) 

Therefore, the Dirac bracket (2.20) contains two separate parts, one due to inverse of 
C^P^ = s-^d the other due to C^J^, defined similarly. So it is enough to do the 

calculations just for the set of $'s. Using the integral form of the constraints: 

/ da6{a)d:^i{a,T) n = 0, 2, • • • 

(2.23) 

= I daS{a)d:^ni{a, r) n = 1, 3, • • • 





n — 


1,3, 




n — 


2,4, 




n = 


1,3, 


I— irr 1 


n = 


2,4, 



it is straightforward to calculate 







n, m 



= 1,3,-- 




mm 



< -2{MB),j J (5((7)(5((7')9™+^9^,(5((7 - a')dada' n, m = 0, 2, • • • 

Mi^- / 5{a)5{a')d^+^d^,5{a - a')dada' n = 1, 3, • • • m = 0, 2, • ■ ■ 



(2.24) 



This matrix possesses a finite dimensional part with i , j indices and an infinite dimen- 
sional part with m, n indices. The finite part can be inverted easily, while for infinite 
part there exist serious difficulties. In fact the most crucial point in studying the prob- 
lem is here. Among so many authors that escaped the existence of infinite number of 
constraints, there are two main references, published almost simultaneously, that have 
tried to invert the matrix C"™" above, and calculate the Dirac brackets. First, the famous 
work of Rcf. [5] in which the authors tried to write down the inverse of C™"' by using 
some undetermined functions i?„m((T', cr") and 5'„j„(cr', cr"), where it is claimed that these 
functions are omitted during calculating the Dirac brackets. However, we think that it 
is not allowed to do calculations for the midpoints a' , o" of the string and thereafter 
consider them at the end points cr = 0, L In other words, the presence of 5{(t) and 5{a') 
in the formula (2.24) is essential. In fact, it may happen that some expressions vanish 
for intermediate points of the string before going to the boundaries. ^ In this way it is 
not clear that in what sense may the functions Rnm{c', c") and Snmic', c") be defined 
appropriately. 

The next reference in this regard is Ref. [9], in which the authors have tried to 
regularize the delta functions in (2.24) and then invert CfJ^. However, in practice this 
method does not seem to make it possible to find directly. The authors have written 
the answer by considering some desired properties of Dirac brackets. In this way the 
main features of the answer is derived, but unfortunately it contains the regularization 
parameter of delta functions, which does not sound plausible. The reality is that, apart 
from some exceptional references [19, 16], most authors have tried to find, in different 
ways, the original results which imply the non-commutativity of coordinate fields 
at the boundaries (see Eqs. (5.23) in the following). We will show in section (5) that 
applying the familiar approach of mode expansion in the context of constrained systems 
gives reliable results. 

Before that we prefer to spend two sections to establish the existence of infinite number 
of constraints by studying the fundamental discrete model corresponding to the continuum 
model. 

^This precise point made the author of [19] to deduce that the first level constraints commute, and 
finally that the coordinate field are commutative at the boundaries after quantization. Being more 
accurate [$(cr, r),$(cr',T)] = for arbitrary a, a' due to antisymmetry of B, while using (2.23) shows 
that at boundary a = 0, [$(0, r), $(0, r)] 7^ as can be seen from (2.24). 



3 Discretization 



The discretized Lagrangian corresponding to the model given in (2.3) can be written as: 

-I N 9 1 N-1 . Y _ V \ 2 N-l 

L- leY. (^n) - ^6 E "^- ) + E X^B - X„) (3.1) 

^ n=0 ^ n=0 V c / 

For the sake of simpUcity in notation we have dropped the i, j indices on 5-field. So B 
should be considered as a matrix and X„ as a column vector in the above equation as 
well as in the following. All associated quantities to Xn carry the same hidden index of a 
column vector. The continuum limit is achieved by the following replacements: 

^ oo 
e ^ 
Ne I 

X^^Xia,T) 

^"+7^" ^9.X(a,r). 

Here wc have ascribed the right difference (divided by e) to spatial derivative. It is also 
possible to do the same with the left difference. In the continuum limit the physical 
quantities in the neighboring points arc the same. 

The Euler-Lagrange equations of motion for intermediate points are: 

Xn-^ + eB^^O 3.3 

where 

An = Xn+l — 2Xn + Xn-1- (3-4) 

The last term in (3.3) is of order e and vanishes in the continuum limit (3.2) giving the 
wave equation (2.5) which implies that the i?-field has no effect in the intermediate points 
and appears only in the equations of motion of the boundary points, as follows 

eXo - - 5(2^0 - Xi) = 

eX^ + ^^H-Asl^ + BX^_, = 0. 
e 

In the continuum limit the terms proportional to e in Eqs. (3.5) vanish, while Xi and 
Xn_i can be replaced by Xq and Xjv respectively, to obtain the following acceleration-free 
equations 

d,X{a,r) + BX{a,r)\,=o,i = 0. (3.6) 

Eqs. (3.6) are the primary Lagrangian constraints. In the Hamiltonian formalism the 
momenta conjugate to coordinates X^ are 

Pn = eX^ + B{Xr,+i - Xn) n = 0, 1, ■ ■ ■ , iV - 1 ^^^^ 
Pn = eX„. 



The canonical Hamiltonian reads: 

1 1 

^ = ^ E [iPn - B{X^+l - Xn)T + {Xn+1 - X„)'] + -p^. (3.8) 
n=0 

To achieve the continuum hmit one should complete the hst given in (3.2) as follows: 
Pn^O n = 0,l,---,N 

^ ^ U{a, r) = X{a, r) + Bd,X{a, r) n = 0, 1, • • • , TV - 1 (3 9) 

^ ^ 7i(n(a, T),X{a, t), d,X{a, r)) 

where in the second line Eq.(3.7) is used and Il{a, r) and H are the momentum field and 
Hamiltonian density, respectively. The canonical equations of motion for the intermediate 
points are as follows: 

p„ = ^[S(p„-p„_i) + MA„] n= l,2,---,iV-l (3.10) 

X„ = ^[pn-^(^n+l-^n)] < n < TV. (3.11) 



For the end points we have 



Po = -JBpo + M{X^-Xo)], 
Pn^-^ [BpN-i + M{Xn - Xn_i)] . 



(3.12) 



The right hand side of the above equations are finite in the continuum limit, while the 
left hand side vanishes (see Eqs.( 3.9)). So Eqs. (3.12) can be viewed as the following 
primary Hamiltonian constraints 

= Md,X(a,T) + BU(a,T)l=o,i- (3.13) 

It should be noticed that (Pn—Pn-i) is of order e^, so in the continuum limit Pm-i/ e can also 
represent the the end point momentum Il{a = I, r) as well as p^/e. For the same reason 
(Xpf — Xiv_i)/e can be interpreted as (9o-X((T, r)|o-=/ despite our previous convention of 
attribution of right difference to the spatial derivative at a given point. These constraints 
can also be derived from Lagrangian constraints (3.5) upon inserting X{a,T) from (2.7). 
The constraints at the end points are completely similar to each other. So without losing 
any point, we can concentrate only on the boundary a = 0. The same arguments can be 
estabhshed for the other boundary a ~ I. 



4 Consistency condition for the constraints 



In the previous section we showed that the equations of motion for the end points of the 
string can be treated as primary constraints. As in any constrained system one should 



investigate the consistency of the constraints. In discrete model this means that one 
should differentiate the equations producing the constraints, i.e. Eqs. (3.12) with respect 
to time to give 

^ [5po + M{X, - Xo)] = 0(e). (4.1) 
Using Eq. (3.11) to insert Xq and Xi into (4.1) gives 

-Bpo + \M{-BA,+pi-po)=0{e). (4.2) 
Inserting MAi from Eq (3.10) in Eq (4.2) results in 

^B{po-pi) + ^{pi-po) = 0{e). (4.3) 

The first term is of order e while the second term is the discrete version of (?o.n((T, t)|o-=o- 
Therefore in the limit e — > the second level constraint emerges as 

$2 = 9,n((7,T)U=o = (4.4) 

One should proceed the consistency condition for the new constraint Due to some 
technical difficulties if we wish to do this in the discrete model it is not a good idea to 
differentiate directly the discrete version ^i^^ of the constraint Instead, it is better 
to represent $^ with In fact we can transfer the condition d(jIl{(7,T)\„=Q = to the 

right by infinitesimal distance e. This is reasonable since in the continuum limit every 
local condition on the fields should be valid in a small neighborhood of a point, not just 
strictly at the given point. This is in fact the "continuity hypothesis". In other words, it 
is not plausible to go to continuum limit just by taking the limits given in equations (3.2) 
and (3.9). It is also needed to impose the continuity hypothesis on the physical quantities, 
which implies that the difference of discrete values of fields in the neighboring points can 
be at most of order e. Therefore, differentiating = + 0(e) gives 

^' = ^(p2-pi) + 0(e). (4.5) 

It is worth noting that differentiation with respect to time does not change the order of 
a quantity. The reason is that time derivative of a quantity is achieved by the Poisson 
bracket of that quantity with the Hamiltonian which is of order e; but in computing the 
Poisson bracket, one differentiates with respect to canonical momenta Pn, which are also 
of order e. So the net result is of the same order. We insert pi and p2 from Eq. (3.10) 
into Eq. ( 4.5), to get 

= ^ [B{p2 - 2pi + po) + M{Xs - 3X2 + 3X1 - Xo)] + 0(e) (4.6) 
Going to the continuum limit we have 



$3 



BdlU{a, r) + MdlX{a, r)] U=o = (4.7) 



where we have used (3.13). 

To investigate the consistency of the strategy is the same as before: we should 
differentiate with respect to time, but we need to transfer the terms one step to the 
right. Then using the equations of motion (3.10), (3.11) and the continuity hypothesis, 
and going back to the continuum hmit, one gets the next constraint as 

$^ = 9^n(a,r)U=o. (4.8) 

It is reasonable that this procedure will produce at the boundaries the infinite set of 
constraints 

$\ a,n, a^n,--- (4.9) 

which is the same as (2.18) and (2.19). One may wonder about the validity of constraints 
not only at the end points, but also in an infinite number of their adjacent points. The 
answer is that a real continues system consists of uncountably infinite points, while in the 
discretized model one imposes the constraints on a countable infinite number of points 
in the vicinity of the boundaries. So physically speaking, nothing bad has happened. In 
other words, suppose we extend the validity of the constraint 9crX(o", r)|o-=o = (in the 
case of a simple string with free end points) to a large countable number of the adjacent 
points of the boundary. Even when the number of points goes to infinity, it will still 
remain in an infinitesimal neighborhood of the boundary. In other words, it will never 
extend in the continuum limit to a finite distance from the boundary. 



5 Reduced phase space 

In this section we try to find out the most suitable basis to describe the physical (re- 
duced) phase space. Whenever second class constraints exist, one should first impose 
the constraints to eliminate the redundant variables and reach the reduced phase space. 
One should then try to find the most suitable bracket on the reduced phase space. It 
is clear that the ordinary canonical quantization procedure (i.e. converting the Poisson 
brackets to commutators) is not consistent in the original phase space, since quantum op- 
erators corresponding to constraints have non-vanishing commutators which contradicts 
the necessity that they should vanish either strongly or on the physical states. However, 
a consistent quantization procedure can be done in the reduced phase space. This is 
achieved by converting the induced brackets on the reduced phase space to commutators. 

Fortunately the famous Darboux theorem ensures us that a unique and well defined 
bracket, which is the same as the Dirac bracket, exists on this space [24]. In fact, the 
Poisson bracket in the original phase space induces the Dirac bracket on the reduced phase 
space [1, 2]. In other words, for any two functions f{q,p) and g{q,p) one can write 

{/(9,P) , 9{(1,p)}d.b = {/(?,P)k=o , 9(q,p)W=o} , (5.1) 

where f{q,p)\<i>=o means evaluation of f{q,p) on the constraint surface described by the 
equations $ = 0. 



In the general case, the second class constraints may be some complicated functions 
of the coordinates. If so, the constrained and physical degree of freedom are mixed with 
each other and it is not generally an easy task to separate them. Sometimes it is almost 
impossible to compute the Dirac brackets directly from the definition (2.20), as is the 
case for our current model (string in the background B-field). Moreover, it may happen 
that the resulted quantum algebra is difficult to handle, specially in order to find the 
corresponding representations. 

Now consider an idealized model in which the second class constraints are given by 
Qk+i, ■ • - Qn , Pk+1, ■ • 'Pn where the coordinates {qi ■ ■ ■ Qn , Pi ■ ■ -Pn) describe the original 
phase space in which the Poisson brackets are defined as 



It is clear that the reduced phase space with coordinates {qi - ■ ■ Qk , Pi - ■ ■ Pk) acquires a 
natural bracket in which summations run from 1 to k, i.e. 



Such a truncated Poisson bracket is a realization of the instruction (5.1) and can be checked 
that is equal to the Dirac bracket of (2.20). In fact, after imposing the constraints on 
this system, the Dirac brackets of the remaining variables are the same as their Poisson 
brackets. Hence, quantization of the system may be achieved in the most simple way, 
such that the familiar algebra oi x — p variables and the corresponding representations 
are still vahd. 

Due to extreme simplicity of the above system, it is much more convenient to change 
the coordinates of a theory with second class constraints to a basis in which the constraints 
constitute a set of conjugate pairs. We call such a coordinate system as normal coordi- 
nates. In most physical theories the Fourier modes are normal coordinates. Traditionally 
people are used to impose assumed commutation relations among the Fourier coefficients 
in order to quantize a field. However, the important point, which is not clearly stated in 
the literature, is that the Fourier modes are the normal coordinates describing the reduced 
phase space. In other words, they carry the "net physics of the theory", hence they are 
independent variables which are suitable for quantization. 

5.1 Free bosonic string 

For example consider a free bosonic string with Neumann boundary conditions [9]. In 
this simple case we are given the conjugate fields X{a,T) and H((j, r) with the following 
Poisson bracket 




(5.2) 




(5.3) 



{X(<j,t),H(<j',t)} = %-<7'). 



(5.4) 



Suppose the fields are real. The most general form of their Fourier expansion can be 
written as: 

1 f°° 

X{a, t) — — / dk [afe(r) cos ka + bkir) sin ka] 

v27r J—oo 

(5.5) 

I roc ^ ^ 

n(o-, t) — / dk [ckir) COS ka + dkir) sin ka]. 

v27r J—oo 

It is obvious that afc(r) and Cfc(r) should be even functions of k while hj.{T) and dfe(r) 
should be odd. Using (5.4) it is easy to see that (0^,0^) and {bk,dk) are conjugate pairs 
in the space of new variables, i.e. 

{akir), Ck'ir)} = {^^(t), 4'(^)} = S(k - k') (5.6) 

and all other Poisson brackets vanish. The second class constraints of the system are: 

9f+iX(<7,r)U=o,. = 

n = 0, 1,2,--- (5.7) 

a2-+in(f7,r)|^=o,/ = 



which can be derived similar to section (2) or by imposing B — Q on its results. The 
constraints at cr = 0, give 

dk{-lfk^''+%k - 

(5.8) 

dA;(-l)"A;2"+^4 = 



00 
00 



/. 
/ 

Since and d^ are odd, Eqs. (5.8) can be satisfied for all n, iff bk — dk = 0, hence 
Eqs. (5.5) change to 

1 

X{a, t) — . — / dkakir) cos ka 
v27r J -00 



n((T, r) = — ;= / dkckir) COS ka . 



This means that in the basis of Fourier modes, b^ and d^ are constrained variables and 
the reduced phase space is simply achieved by omitting them. This is the advantage of 
using the Fourier modes as normal coordinates. If we were insisting on working in the 
original basis X{a, r) and II (a, r), we would encounter difficulties explained at the end of 
section (2). Now let us consider the constraints on the end point a — I. They give 

/oo 
ciA;(-l)V+^afcsinM = 

(5.10) 

dA;(-l)"it2"+^CfcSinM = 



00 
00 



Since and arc even, the integrands in Eqs. (5.10) are even with respect to k. So the 
only way to impose the constraints is: 

ak — Ck — ^ for /cZ 7^ 0, TT, 27r, ... (5-11) 



Once again we sec the miracle of working with Fourier modes. In this basis a large class 
of the coordinates and are omitted due to the constraints, just remaining with those 



with discrete values ior k as k — nn/l for integer n. Finally the original field variables 
can be expanded in terms of infinite countable Fourier modes a„ and Cn as canonical 
coordinates of the reduced phase space as: 

(5.12) 

V! V T ^ ~r 

It is easy to see that the Fourier modes am and obey the canonical brackets: 

{(^mi (^n} — {Cm; C„} = 0, {Orn; C„} = (^mn (5.13) 

Using expansion (5.12) and brackets (5.13), one can compute the Dirac brackets of any 
two physical functions of the original variables X(o", r) and n(o", r). This means that we 
have followed the prescription given in Eq. (5.5) to find the Dirac brackets. As is well 
known [1], the second class constraints should strongly vanish before quantization. This 
fact can be stated in terms of the Fourier modes more clearly. The constrained modes 
{bk,dk) for all k and {ak,Ck) for k ^ nn/l should vanish before quantization. Then one 
can quantize the theory by assuming canonical commutation relations among {am, Cm) in 
the expansion (5.12) as 

[om, c„] = iM^n- (5.14) 



5.2 Open string in background 5-field 

Let us now consider the string in the background 5-field. We expand the main fields 
X(cr, r) and n((7, r) as done in (5.5) with (afc,Cfc) and {bk,dk) as conjugate pairs. The 
constraints (2.18) and (2.19) can be rewritten as 

af$(a,r)U=o,/ = 

n = 0, 1,2, ••• 5.15 

af+^n(a,r)U=o,. = o 

Using the constraints (5.15) at cr = read 

/oo 
dk{-lYk'^''{kMhk + Bck) = 

(5.16) 

dA;(-l)"A;'"+'(4) = 0. 



oo 
oo 



Remembering that (6^, dk) are odd and Ck is even with respect to k, the constraints (5.16) 
are satisfied for all n and A; 7^ iff 

bk = —M-'Bck 

k (5.17) 

4 = 0. 



Now imposing the constraints (5.15) on the end point a = I gives 

dk{-Vfk^''{-kMa^,) sin kl = 



/ 

;r (5.18) 

/ rffc(-l)"fc'"+'(cfc)sinA;/ = 0. 

J—oo 

These equations show that {ak,Ck) and consequently bk should vanish for k ^ mr/l {n 
integer). To this end, some care is needed for the zero mode {k — 0). In the limit /c — > 0, 
using (5.17) we have 

lim sin /ccr = lim ( -jM'^Bck sin ka) = -M~^Bcoa. (5.19) 

Therefore the linear term {—M^^BcQa), coming from the zero mode of the sine terms 
should be present in the expansion of X{a,T). Similarly to (5.12) the term is also 
necessary as the zero mode of cosine terms. However according to the global symmetry 
given in (2.4), we are allowed to add any constant term to the expansion of X{a, r). This 
term should not disturb the validity of constraints (5.15) and should vanish in the limit 
B ^ 0. We fix this arbitrariness by adding the constant term as M^^BcqI/2. As we will 
see later, Cq is constant according to the equations of motion and moreover, this choice 
guarantees that the coordinates of the center of mass of the string are commutative. 

Putting all these results together, the most general form of the fields satisfying the 
constraints can be written as 

X(a,T)^^ [ao-M ^Bco{a - -) \ + d - >J a„cos— M ^Sc„sm- 



VI \ 2 y V ^ ^1 V I mr I 



(5.20) 

These relations show that in the case of mixed boundary conditions again a„ and c„ are 
suitable canonical coordinates of the reduced phase space. Note that (a„, Cn), as canonical 
coordinates, still obey the canonical brackets (5.13). Hence, from the general prescription 
given in (5.1) it is easy to calculate the Dirac brackets of the fields, just by using their 
expressions in terms of normal coordinates a„ and c„. As can be seen from (5.20) the 
S-field appears only in the expansion of coordinate fields Xi{a,T), while the momentum 
fields Hj((T, r) are unchanged. Therefore, the Dirac brackets {Xi{a,T),Ilj{a' ,r)}D.B and 
{Hi((7, r), Hj((7', t)}d.b are the same as the corresponding Poisson brackets given in Eqs. 
(2.13). However, for the Dirac brackets of coordinate fields, from (5.13) and (5.20), one 
finds 



{X,{a,T),Xj{a',T)}D.B = {M-'B), 



1 + - 2^ -sm— fa + a" 



I T^^in I 

The summation over the sines is the Fourier expansion of saw waves as follows 



(5.21) 



oo 



y -sinne ^ I ~ 5.22 

] Un-e) 0<9<7r 



2^ 



This function is discontinuous at ^ = 0, 27r, .... Supposing its values at these points to be 
the average of right and left limits, i. e. zero, we can write the final result as: 

{Xi{a, T),Xj{a', t)}d.b = a,a'^0 
{Xi{0,T),Xj{0,T)}D.B = -2{M-'B)ij, 
{XS: r),Xj{l, r)}D.B = 2{M-^B)ij. 

As we see, after quantization the coordinate fields are noncommutative at the end-points 
of the string, in agreement with the well-known results given in the literature [4, 5, 9]. If 
we had not added the constant term M~^BcqI/2 to the expansion (5.20) the above result 
would have differed from (5.23) just by a constant term throughout the string, as well as 
at the end-points. In other words the non-commutativity at the end-points (7 = and 
(7 = I have opposite signs since we have imposed the condition that the center of mass 
coordinates are commutative. 

Our emphasize in our derivation of the important result (5.23) is that we have not 
used the expansions of fields in terms of the solutions of the equations of motion. In fact, 
we have not considered the time dependence of the physical modes a„(T) and c„(r) which 
should be determined by means of the special form of the Lagrangian or Hamiltonian. 
This feature of our approach will be explained more in the next section. 

6 Quantization and equations of motion 

The traditional canonical quantization procedure is as follows: one considers the general 
solution of the classical equations of motion, then imposes the boundary conditions to de- 
crease the number of possible modes ^, and finally assumes suitable commutation relations 
amongst the physical modes to quantize the theory. One may ask: "is it really necessary, 
or even allowed, to use classical equations of motion in the process of quantization?" . 

Let us first study the problem in an ordinary (unconstrained) system. The important 
point is that the special form of the Hamiltonian (or Lagrangian) is not the essential point 
that determines the algebra of physical observables and consequently the structure of the 
space of physical states of the theory. On the other hand, the role of the Hamiltonian 
is just determining the dynamics of the system. Given the initial state of a system, the 
Hamiltonian is the main tool which gives the time evolution of the state of the system. 
However, people usually try to construct the basis of the space of the physical states 
in such a way that the Hamiltonian operator is diagonal; since this provides consequent 
convenience to follow the dynamics of the system. 

For example in one dimension whatever the Hamiltonian is, one can use the algebra of 
the X — p operators acting on the space with \x' >'s or \p' >'s as the basis, as well as the 
equivalent algebra of a — operators and the corresponding basis in which the operator 

*The classical equations of motion are usually linear differential equations; so one can expand their 
solutions using a complete set of modes 



(5.23) 



a'a is diagonal. However, it is well-known that for a free particle the former algebra is 
more suitable while for the harmonic oscillator the latter is more appropriate to study the 
dynamics of the system. 

To this end, we want to emphasize that conceptually it is not needed to treat the 
quantum fields as the solutions of the classical equations of motion. However, the reader 
may have encountered several books or papers where the authors write down the fields 
as an expansion in terms of solutions of the classical equations of motion (for example 
plane waves with definite uj — k relations) and then quantize the theory by imposing 
assumed canonical algebra among the coefficients of the expansion. A careful notice leads 
to the observation that most of the time, the explicit time dependence of the terms in 
the expansion are not used during the subsequent analysis. For example in harmonic 
oscillator problem, it is just the quantum algebra [a, a^] = 1, coming from the classical 
algebra {x,p}pB. = 1 upon quantization, which determines the basis of physical states as 
|n>,|n+l >,■■■. ^ Then the explicit form of the Hamiltonian may be used to determine 
the time dependence of a{t) and o){t) in Heisenberg picture as a{t) = a(0) exp{—iujt) and 
a^(t) = (0) exp[iut) , which is not essential in determining the quantum properties of 
the observables as well as the space of physical states. 

Our experience shows that in quantum field theory the Fourier expansion of the fields 
come out, most of the time, to be useful in the process of quantization. As discussed 
in the case of the models considered in this paper, this is just a change of variables in 
the phase space from X((T, r) and H(ct, r) say, to ak{T) , bk{T) , ■ ■ ■ etc. Then, regardless 
of the dynamics of the system, the Poisson brackets of the original variables determine 
those of the new ones. There are two main advantages in this change of variables. First, 
the Hamiltonian may be diagonal or have a simpler form in the new basis. Second, the 
constraints as well as boundary conditions (which are also considered as constraints in 
our approach) may be applied in a simpler way in the framework of the new variables. 
Therefore, the Hamiltonian has some partial role in quantization since the dynamics 
of constraints, (i.e. the consistency of the constraints) should be investigated classically 
before quantization. It is not possible to construct a quantum algebra among the variables, 
unless the Hamiltonian has vanishing brackets with the constraints on the physical space 
(reduced phase space). ^ 

In other words, although the full dynamics of the physical variables is not essential 
for quantizing the theory, the dynamics of the constrained variables should be worked 

^Note that imposing the condition of unitarity on the physical states, restricts n to positive integer 
values, which in this case guarantees that the energy states are bounded from bellow. Therefore, besides 
the quantum algebra of observables, some other physical requirements such as unitarity principle play 

important roles in determining the set of physical states. 

^If the constrained quantities are assumed as vanishing operators in the quantization procedure, then 
their brackets with the Hamiltonian should also vanish. On the other hand, if one quantizes the theory 
by imposing the condition that quantum operators corresponding to constraints should kill the physical 
states, then again it is clear that their brackets with the Hamiltonian should also kill the physical states. 



out completely before quantization, so that the final brackets of the Hamiltonian with 
the constraints vanish. This means that using the classical equations of motion, all the 
secondary constraints should be computed before quantization. Schematically we can say 

primary constraints + classical equations of motion — > secondary constraints. (6.1) 

This point can be seen clearly in the example of the bosonic string with Neumann bound- 
ary conditions. If one had considered just the primary constraints daX{a,T)\a=o,i — 
instead of the whole set of (5.7), then one would not have been able to deduce the expan- 
sion (5.12) for the fields. However, if one considers the primary constraints together with 
the equations of motion (resulting from the Hamiltonian (2.8) with B — 0) one obtains 



drX{a,r)=U{a,T) 
drU{a,T)=dlX{a,r). 

Then one can easily check that 

drid^Xia, r)) = d^idrXia, r)) = d^i^, r) 

In this way the infinite set of constraints(5.7) are in fact resulting from the primary 
constraints plus the equations of motion (see (6.1)). This argument shows that there is 
no way to escape the fact that an infinite number of constraints really exist. If one uses 
the full capacity of the classical equations of motion (6.2), or the explicit form of the 
Hamiltonian in terms of the normal coordinates as 

^ TT/TT 

H ^ -J2^nCn + 7;J2^riani-ry, (6.4) 
^ n=0 ^ n=l ' 

in order to determine the dynamics of the physical variables an(T) and c„(t), then one 
obtains d„ = c„ and c„ = — (n7r/i)^a„, which acquire the solution 

/ nn \ I / nn \ 

an{r) = a„(0) cos — r H c„(0) sm — r 

(6.5) 

Cn(T) c„(0) COS {'Y'T^ - -J- «n(0) sin {-y'T^ ■ 

However, we insist again that the full dynamics of the physical variables is not necessary 
to quantize the theory. It seems that this partial role of the classical equations of motion 
in determining the dynamics of the constraints (or boundary conditions in most of the 
familiar physical systems) is the hidden reason behind the common practice of expanding 
the fields in terms of the classical solutions of equations of motion before quantization. 

Now let us consider again the string in background i?-field to observe the above points. 
In this case the equations of motion resulting from the Hamiltonian(2.8) read 



drX{a, t) = H((7, r) - BdaX{a, r) 
drIl{a,T) = Bd^U{a,T) +Md^X{a,T). 



(6.6) 



This gives 



drma,T)) = d„{U{a,T)) 
dridM^,r))=dima,r)) 



(6.7) 



In this way the fuU set of infinite constraints (5.15) emerge as the resuh of combining 
the primary constraints $(0", r)|o-=o,/ = with the equations of motion (6.6). Using the 
expansion (5.20), the canonical Hamiltonian (2.8) in terms of the normal coordinates can 
be calculated as 

H^-Y. c^M-'c^ +-Y: anMa^C^f. (6.8) 

^ n=0 ^ n=l ' 

We observe again that the full content of the dynamics is not needed for the quantiza- 
tion process. In fact, using the Hamiltonian (6.8), or the equations of motion (6.6), to 
determine the dynamics of the physical variables a„(T) and c„(t), one finds 

d„(r) = M-ic„(r) 

/nvTN^ n = 0,l,--- (6.9) 



Cn(T) = -(^— ) Ma„(T) 



which acquire the solution 

f rni \ / , ,_i { nn 

a„(r) = a„(0)cos — r H M c„(0)sm — -) 

\ I y nn \ I 



(6.10) 



/ nn \ nn / nn 

Cn{r) = Cn(0)cos ( — T I ^Ma„(0)sm ( — t 

and 

ao(r) = M-ico(0)T + ao(0) 

(6.11) 

Co(r) = Co(0). 

As is apparent in order to quantize the theory, specially finding the important results of 
(5.23), one does not need to know the explicit time dependence given in Eq. (6.10). Also 
note that cq is constant, in agreement with our previous trick of adding the constant term 
M'^BcqI/2 to the expansion of coordinate field in Eq. (5.20). 



7 Conclusion 

In this paper we discussed different aspects of the idea of considering the boundary condi- 
tions as Dirac constraints. Our theoretical laboratory for this aim was " an open string in 
a background i?-field". We observed that besides the singularity of the Lagrangian, the 
boundary conditions can serve as a source of introducing the primary constraints. Ana- 
lyzing in detail the discretized version of the model shows that the primary constraints are 
the continuum limit of the equations of motion for the end points; while the secondary 
constraints are derived by imposing the consistency conditions on these equations and 



then going to the continuum hmit. In this process the continuity hypotheses is impor- 
tant. This imphes that in order to get a continuum solution, the fields in the adjacent 
points in the corresponding discrete model should not differ drastically. In this way it 
turns out that the discretized model highly supports the existence of infinite chains of 
second class constraints. 

The continuity hypotheses is also deeply related to Fourier expansion, as follows. It is 
well-known that in writing any field as a summation over the set of well-behaving contin- 
uous sine and cosine functions, any finite discontinuity in the field or in a finite number 
of its derivatives, is not seen by the expansion and is somehow removed from the prob- 
lem. However, the boundary conditions should not be considered as such discontinuities. 
Although boundary conditions imposes some restrictions just on definite points at the bor- 
der of the medium (i.e. end points of the string), they have their considerable influence 
throughout the whole medium. The important point is that the Fourier expansion plays 
the role of a carrier of boundary conditions from the boundaries through the medium. In 
fact, the emergence of an infinite number of constraints causes serious restrictions on the 
Fourier modes invited to the expansion of the fields. Since the Fourier modes are alive 
in the whole medium as continuous and well-behaved function, the message of boundary 
conditions is distributed in this way throughout the system. The familiar example in 
this regard is an open string with Neumann boundary conditions, in which the fields are 
expanded in terms of a set of discrete cosine modes only. In this case the Fourier expan- 
sion is, in fact, used to soften and fiatten the fields undergoing definite conditions at the 
boundaries. 

Summarizing, the continuity hypotheses implies, in the discrete model, the validity of 
constraints should be spread in a set of countable infinite number of adjacent points near 
the boundaries. This fact then shows itself in emerging infinite number of constraints; and 
finally causes omitting a large class of Fourier modes, which leads somehow to propagation 
of the effect of boundary conditions through the medium. 

Another aim of this paper was studying the Poisson structure of the reduced phase 
space. Using the original coordinates of phase space implies serious difficulties in com- 
puting the complicated Dirac brackets. We observed that Fourier modes can serve as the 
normal coordinates of the reduced phase space. Using the Fourier modes make it possible 
to do calculations. Moreover, it gives a valuable understanding of the Poisson structure 
of the reduced phase space which contains the true physical degrees of freedom of the 
model. In this way after disappearing some Fourier modes as redundant (constrained) 
coordinates, the remaining modes can be viewed as the physical degrees of freedom, i.e. 
the coordinates of the reduced phase space. Fortunately these modes emerge as conjugate 
pairs with a simple and well-established bracket. As explained in the text, the brackets 
of the remaining modes define the Dirac brackets of fields. In this way one expands all 
associated fields and quantities just in terms of the normal coordinates, and then using 
their brackets one writes down all the Dirac brackets. It turns out that for an open string 
with mixed boundary conditions there remains no doubt about the fact that the Dirac 



brackets of the coordinate fields at the boundaries of the string are nonzero due to the 
S-field. Then, upon quantization, the coordinates of the string, and hence the coordinates 
of the D-brane, are non-commutative. 

Another new feature in our approach is that we do not use of solutions of equations 
of motion in the process of quantization. We showed clearly, in the model under consid- 
eration, that it is in fact possible to find the algebraic structure of the quantum theory 
without any need to expand the fields in terms of solutions of classical equations of mo- 
tion. We argued that in any quantum theory, constructed upon quantization of a classical 
model, one only needs to consider the dynamics of the constraints before quantization. 
In other words, it is not essential to find the dynamics of the complete set of physical 
variables to do quantization. 

We think that our approach may be useful in analyzing the physical structure together 
with quantization of every model with complicated boundary conditions. Two recent 
examples can be seen in [25, 26]. This approach may be applied also to more complicated 
systems such as membranes [27, 28, 29]. 
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